Spin Torques in Ferromagnetic/Normal Metal Structures 
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Recent theories of spin-current-induced magnetization reversal are formulated in terms of a spin- 
mixing conductance G""". We evaluate G™*^ from first-principles for a number of (dis)ordered 
interfaces between magnetic and non-magnetic materials. In multi-terminal devices, the magnetiza- 
tion direction of a one side of a tunnel junction or a ferromagnetic insulator can ideally be switched 
with negligible charge current dissipation. 
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"Giant magnetoresistance" refers to the large change 
of resistance brought about by applying an external mag- 
netic field so as to change the angle between the magne- 
tization directions of adjacent magnetic films H. Since 
a spin injected into a magnetic material experiences a 
torque, it has been argued that there should be a re- 
verse effect, namely, that passage of a current through 
adjacent magnetic layers should lead to the transfer of 
spin angular momentum from one layer to the other [yj^] 
with possible reorientation of the magnetizations for suf- 
ficiently large currents [g| . Interestingly, the sign of the 
corresponding torque should be reversed on changing the 
current direction leading to the possibility of making 
an electronically accessible non-volatile magnetic mem- 
ory whose performance on downscaling would appear to 
compare favorably with other alternatives Q. The ex- 
perimental observation of current-induced magnetization 
reversal ( "spin transfer" ) effects by a number of groups 
|g|-|[ calls for a more general framework for electrical 
transport in non-coUinear and disordered magnetic sys- 
tems. Two suitable and nearly equivalent approaches are 
the circuit theory based on kinetic equations for spin cur- 
rents [ pJ3[ and the random matrix theory of transport in 
magnetic heterostructures ]ll[ . 

In the circuit theory, the torque and current are formu- 
lated in terms of (real) spin-up and spin-down conduc- 
tances G^ and G^ and a new spin-mixing conductance 
Qmix •^\^\q\^^ jg complex. The mixing conductance has 
only been studied using simple free electron models and 
little is known about its dependence on real material pa- 
rameters. Since free-electron models miss an important 
contribution to spin transport in layered magnetic mate- 
rials coming from the mismatch of the complex d bands 
responsible for itinerant ferromagnetism [|l2| , [13| , it is im- 
portant to take this into account when evaluating G™" 



for realistic interfaces. 

In this paper we use methods recently developed 
|l#-M to calculate the scattering matrix within the 
framework of density functional theory in order to evalu- 
ate the spin- and spin-mixing conductances for a number 
of systems of current interest. We show that the spin- 
mixing conductance in tunnel junctions can remain large 
even when the conventional conductance itself is made 
vanishingly small, which has not been anticipated pre- 
viously and should be important in utilizing the spin- 
torque effects. We suggest a configuration in which this 
effect could be observed. 

We begin with two basic elements of a magnetic circuit: 
a non-magnetic (NM) element in which there is a spin- 
accumulation in the direction s [as a result of injection 
from elements of the circuit e.g. by a source-drain current 
in Fig. 11(a)] and a ferromagnetic (FM) element whose 
magnetization is given by the unit vector m. In general 
s and m are non-collinear and it is convenient to split the 
2x2 matrix current / from the non-magnetic element into 
the ferromagnetic element into a scalar charge current /q 
and a vector spin current Igi / = (/o -f <t • Is)/2 where a 
is a vector of Pauli spin matrices. 



I, = m[(GT - G^){f^ f^) + [G^ + G^)^f 
+ {G^ + G^ - 2ReG"")s • mA/^] 
-2sReG™"A/^ + 2(s x m)ImG™"A/^ 



(1) 



and J^ ^ J^ +a- sA/^ and /^ = f^ + ct ■ mA/^ 
describe the deviation of the distribution functions in 
the NM and FM nodes, respectively, from their equilib- 
rium values p7[ . When a spin-current is injected into a 
ferromagnetic material, the component of I^ perpendic- 
ular to the magnetization direction m (times the Bohr 
magneton and gyromagnetic ratio) equals the torque act- 



ing on the ferromagnet ||ll|] and this is seen to be de- 
termined entirely by the real and imaginary parts of 
Q-mix ^ The spin-up, spin-down and spin-mixing con- 
ductances characterizing transport through an interface 
are defined as G^ = ^ir(ijt|), G^ = ^tr{t\ti) and 
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transmission and reflection matrices t 



rlr|), in terms of the spin-dependent 
and r-^ny, I is 

an AI X M unit matrix where M is the number of conduct- 
ing channels in the NM element |l^Jll|]. Note that the 
mixing conductance can be measured in principle by the 
angular magnetoresistance of perpendicular spin valves 
M analyzed by circuit theory |19|| . 

The parameter- free calculation of the transmission and 
reflection matricesJM is based on the surface Green's 
function method |20[ implemented with a tight-binding 
linear muffin tin orbital basis pit]. The calculation time 
scales linearly with the number of layers in the scatter- 
ing region and as the cube of the size of the in-plane 
unit cell. Because a minimal basis set is used, we are 
able to perform calculations for lateral supercells con- 
taining as many as 10 x 10 atoms and to model disorder 
very flexibly within such supercells without using any 
adjustable parameters. The electronic structure is de- 
termined self-consistently within the local spin density 
approximation. For disordered layers the potentials are 
determined self-consistently using the layer CPA approx- 
imation |20| . The calculations are carried out with a kn 
mesh density equivalent to more than 3600 mesh points 
in the two-dimensional Brillouin zone (BZ) of a 1 x 1 in- 
terface unit cell. The results of calculations for clean and 
disordered fcc(lll) Cu/Co and bcc(OOl) Cr/Fe metallic 
interfaces and for an fcc(lll) Cu/Co/Vac/Co tunneling 
configuration are given in Table 0. 

We first discuss our results for the Cu/Co and Cr/Fe 
interfaces. Both have been the subject of much study in 
the context of exchange coupling and giant magnetore- 
sistance and different calculations of the interface trans- 
mission matrices yield very similar results for the spin- 
dependent interface resistances |l4|-|l6| . The atomic vol- 
ume of each pair of materials is very similar and we model 
disordered interfaces using a 2 layer thick 50-50 alloy of 
the component materials. The results given in the Ta- 
ble do not depend sensitively on the alloy concentration 
used. 

For a clean Cu/Co(lll) interface, the real part of the 
mixing conductance is comparable in size to the spin-up 
and spin-down conductances but the imaginary part is 
almost a factor of 50 smaller. Interface disorder increases 
the mixing conductance, the real part by about 35%, the 
imaginary part by a factor of three. 

It is interesting to compare IG™'^] with G^ + GK The 
mixing conductance mainly contributes to the torque 
while G^ +G^ determines the electron current. Large val- 
ues of |G""^|/(G1^ +G^) mean more torque per unit cur- 
rent. The calculations show that disorder at the Co/Cu 



interface increases the spin-torque. Another possibility 
to increase the ratio of jG™"! to G^ -I- G^ is to insert an 
impurity layer on the Cu side; a Co monolayer inserted 
on the Cu side scarcely changes the mixing conductance 
but reduces the normal conductance significantly. 

For the Cr/Fe interface the band structure matching 
and the effect of interface disorder are quite different 
compared to Cu/Co. Whereas the majority spin states 
of Co/Cu match very well, it is the minority spin-states 
in Cr and Fe which match best; see Figs. H(c,d). For a 
perfect Cr/Fe(001) interface, the mixing conductance is 
almost twice as large as the normal conductance. The 
expression for the mixing conductance shows that hav- 
ing a large number of propagating channels on the non- 
ferromagnetic side of the interface can lead to a large 
mixing conductance. 

It is of interest to have a closer look at the k|| resolved 
mixing conductance G™" . We can see from Fig. |l|(c) 
that close to the center of the Brillouin zone (BZ) the 
real part of the G™" is very large, even larger than the 
number of channels in Cr at the same k||-points, shown 
in Fig. 0(b); at the same kii -points the transmission of 
majority spin electrons is very low. Thus at some kii the 
mixing conductance can be much larger than the normal 
conductance. This can be understood in terms of a sim- 
ple one-dimensional infinite barrier model in which the 
spin-up and spin-down barriers are displaced in space by 
an amount A. For both spins the reflection amplitude is 
1. However, the displacement A introduces a phase shift 
g-2ifcA £qj. electrons with wave-number k so that the mix- 
ing conductance can have any value Go(l ^ g"^*'^'^) for 
this simple single-channel model. Although G™" is large 
around kn = for Fe/Cr, the minority-spin reflection is 
very low in most parts of the BZ so that after averaging 
over the BZ, G™*^ is not very high compared with the 
normal conductance. 

The imaginary part of G™" is related to the spin pre- 
cession which results from the non-coUinear alignment of 
the spins of the injected electrons and the magnetiza- 
tion (or an external magnetic field). A non- vanishing 
imaginary part of the mixing conductance, ImG™*^, 
should result in antisymmetry with respect to time re- 
versal |2^j2^. However, ImG™*^ is small in all the sys- 
tems we have studied. The reason for this can be un- 
derstood by examining the kn resolved imaginary part 
of G™'"^ shown in Fig. |l|(f). ImG™" can be negative 
as well as positive and it is the partial cancellation of 
these contributions which leads to the net result being 
small. This can be illustrated using the simple phase- 
shift model as follows. Suppose that the phase shift 
of the reflected waves (5(k||) is distributed randomly be- 
tween ifi and ip2 with equal weights and that the ampli- 
tude of rlri is A, then the average mixing conductance 
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-(p2 thenlniG""^ is zero and 



Re G™" — Gq. In a more realistic treatment the weights 
will not be homogeneous and Lpi will be slightly different 
from —ip2 so that ImG™'^ is not exactly zero. To obtain 
a large ImG™'^ experimentally, one should use a mate- 
rial with a small number of conducting channels such as 
a semiconductor as the non-magnetic element. Also in a 
ferromagnetic insulator Im G™" could be significant psj . 

We may conclude that a large reflection amplitude does 
not mean that the mixing conductance must be small. 
It provides us with a means to realize large values of 
1G™'^|/(GT _I_ Qi-\^ j-,y ygjjjg g^ tunnel junction or ferro- 
magnetic insulator as the FM element. For such non- 
conducting interfaces we predict that it will be possible 
to obtain a non-zero spin current while the electronic 
current is zero or vanishingly small. We confirm this by 
calculating the mixing conductance for a Cu/Co/Vac/Cu 
tunnel junction, here we use 6 layers empty spheres to 
model the barrier. The mixing conductance is found to 
be 0.41 - i 0.04 and 0.53 - i 0.003 for clean and dirty sys- 
tems, respectively, in units of 10^'^J7~^to~^, which very 
close to the mixing conductances of Cu/Co interfaces. 
The normal conductance is of the order of 10~^^ in the 
same system. For even thicker layers the torque simply 
equals that of the metallic interface and the tunnel junc- 
tion only suppresses the charge current. 

To inject a spin-current in the absence of an elec- 
tron current we need to consider a three terminal device 
("spin-flip transistor") such as that sketched in Fig. H(a) 
[|^,|^. A current from FMl into FM2 induces a spin- 
accumulation in the NM node which is easily calculated 
p4| . Here we simply assume that a spin-accumulation 
exists in NM and analyze what happens when FM is a 
magnetic insulator or is the top magnetic element of a 
magnetic tunnel junction. In the latter case it is pos- 
sible to independently determine the orientation of FM 
by means of the tunneling magnetoresistance ||2^ . The 
spin torque is that of the metallic junction, but without 
the energy dissipation caused by the particle current. In 
practical memory devices it may be advantageous to be 
able to achieve this separation of particle and spin injec- 
tion. 

In summary, we have studied the mixing conductance 
of Cu/Co, Cr/Fe and Cu/Co/Vacuum/Co configurations 
taking the full transition metal electronic structure into 
account and including disorder. The effect of the mixing 
and normal conductances can be separated for a three 
terminal device where one of the terminals is a ferromag- 
netic insulator or a magnetic tunnel junction for which 
the normal conductance can be made vanishingly small 
without affecting the size of the mixing conductance. 
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FIG. 1. (a) Sketch of a three-terminal device where a nor- 
mal metal (NM) element is connected to three ferromagnetic 
elements FMl, FM2 and FM. An applied bias causes a current 
to flow between FMl and FM2. If FM is a tunnel junction 
or a magnetic insulator then the particle flow into FM should 
be vanishingly small, (b) number of propagating channels in 
first Brillouin Zone for bulk Cr. For a clean Cr/Fe interface 
(c) and (d) show the minority-spin and majority-spin conduc- 
tances, (e) and (f ) the real and imaginary parts of the mixing 
conductance. Units of conductance are e^ /h. The result of 
integrating over the whole Brillouin Zone is given in brackets 
at the top of each panel. Note that the values represented by 
color differ per panel but with the zeros are represented by 
the green color. 
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